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Abstract 

We apply the Sklyanin method of separation of variables to the reflection algebra 
underlying the open spin-| XXX chain with non-diagonal boundary fields. The 
spectral problem can be formulated in terms of a TQ-equation which leads to 
the known Bethe equations for boundary parameters satisfying a constraint. For 
generic boundary parameters we study the asymptotic behaviour of the solutions 
of the TQ-equation. 
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1 Introduction 

Exact solutions of spin chain models have provided many insights into the properties of in- 
teracting many-body systems subject to strong quantum fluctuations. Various methods have 
been established to study the spectrum and nature of their low-lying excitations of such mod- 
els as well as their thermodynamical properties and even correlation functions without the 
need to revert to perturbational approaches. At the same time various problems concerning 
systems with open boundaries are still not solved completely. Even for the prototype spin-^ 
XXZ chain with general open boundary conditions techniques for the solution of the spectral 
problem have been developed only recently [1,4,11-14]. This model, apart from being the 
simplest starting point for studies of boundary effects in a correlated system, allows to inves- 
tigate the approach to a stationary state in one-dimensional diffusion problems for hard-core 
particles [7, 8] and transport through one-dimensional quantum systems [5] . Its Hamiltonian 
is given by 

L-l 

H = ^2 [ a j + a j°j+i + (°j ^J+i ~ cosn V ~ cosn V 

+ cxf coth £~ H — — (of cosh 6~ + icrF sinh 9~) sinh n 

I smh£- v J 

2k+ 



+ 



a z T coth £ + H (cr? cosh 6 + + \a v T sinh 9 + ) sinh ri 

sinh v ' J 



where a", a = x,y,z denote the usual Pauli matrices at site j and the parameters k ± , 
8 ± refer to the left and right boundary respectively. The Hamiltonian fll.ip is a member of a 
commuting family of operators generated by a transfer matrix based on a representation of 
Sklyanin's reflection algebra [16]. While this establishes the integrability the actual solution 
of the spectral problem through application of Bethe ansatz methods has been impeded by 
the absence of a reference state, such as the ferromagnetically polarized state with all spins up 
for the case of diagonal boundary fields. Interestingly, by imposing certain constraints obeyed 
by the left and right boundary fields, the eigenvalues of the spin-| XXZ chain and of the 
isotropic spin-,!? model can be obtained by means of the algebraic Bethe ansatz [4, 12]. In an 
alternative approach, Nepomechie et al. have been able to derive Bethe type equations whose 
roots parametrize the eigenvalues of the Hamiltonian for special values of the anisotropy 
f] = vk/(p + 1) with p a positive integer and where the transfer matrix obeys functional 
equations of finite order [13]. Their approach relies on the periodicity of the underlying 
trigonometric i?-matrix of the model which is missing in the rational limit rj — > of the 
isotropic chain. For generic values of the anisotropy the spectral problem has been formulated 
as a TQ-equation assuming that the large j-limit of the transfer matrices with spin j in 
auxiliary space exists [19]. No such constraints are needed in the derivation of a different set 
of recursion relations for the diagonalization of fll.ip based on the representation theory of 
the g-Onsager algebra [1]. Very recently, Galleas has formulated another functional approach 
to determine the eigenvalues of fll.ip in the generic case without a reference state [11]. Taking 
certain matrix elements of the transfer matrix involving the ferromagnetic pseudo vacuum 
and the unknown eigenstates he derives the Bethe equations of Nepomechie et al. without 
the need to restrict the anisotropy to roots of unity. Unfortunately, this approach yields no 
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information on the eigenstates, therefore further studies of correlation functions are out of 
reach for now. 

In this paper we approach the problem by means of a different method which circumvents 
the difficulties of the algebraic Bethe ansatz in the absence of a reference state. Within 
Sklyanin's functional Bethe ansatz (or separation of variables method) [17] the eigenvalue 
problem is formulated using a suitably chosen representation of the underlying Yang-Baxter 
algebra on a space of certain functions. This approach was independently confirmed by use 
of Manin matrices [6] and has proven its strength in particular in models with non-compact 
target space lacking a reference state, e.g. the Toda chain [15] or the sinh-Gordon model [3,18] 
where the spectral problem can be formulated in terms of separable functional equations on 
this space which then have to be solved by exploiting the analytical properties of the problem. 

The article is organized as follows. In sections [2] and [3] we will review the basic properties 
of the reflection algebra and its representation for the open XXZ spin chain. Then in section 
H] we will implement the functional Bethe ansatz for the restriction to the open XXX spin 
chain following Sklyanin's original article [17]. Finally in section [5] we address the problem 
of extracting the solution of the spectral problem from the resulting second order difference 
equations and study the analytical properties of their eigenfunctions. The last section is 
devoted to a concluding summary. 



2 Integrable Boundary Conditions 

Sklyanin's construction [16] of integrable systems involving boundaries is valid for a gen- 
eral class of integrable systems characterized by an i?-matrix of difference form R(X,p) = 
R(X — A*)e End(V <8> V) (V is a vector space with dimV^ G N ) which not only satisfies the 
Yang-Baxter equation 

R 12 {\ - p) R 13 (\ - v) R 23 {p -v)= R 23 {p - v) R 13 {X - v) R 12 {\ - p) (2.1) 

but also several conditions such as symmetry with respect to the permutation operator P on 

V <g> V (Px® y = y<g)x), 

R(X) = PR{X)P , (2.2) 
unitarity involving some complex function p(X), 

R(X)R(-X) = p{X) (2.3) 

and crossing unitarity for another complex function p(A), 

R h {X)R h (-A - 2rj) = p(X) . (2.4) 

The parameter r\ characterizes the i?-matrix and the superscript tj denotes the transposition 
with respect to the jth space in the tensor product V ®V . Here we will need the well-known 
6- vertex model solution 



R{X,p) 



A 0\ x _ sinh(A - p) 

b(X,p) c(X,p) 1 

c(X,p) b(X,p) 

\0 1/ 



sinh(A - p + v) , 2 ^ 
sinh?7 

c ( A >/ i ) = ■ ryr — -r 

sinn (A — p + 7]) 
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of the Yang-Baxter equation (12. ip . It generates the Hamiltonian of the spin-| XXZ chain 
with 

, , A i sinh A sinh(A + 2ri) 

= 1 ■ P(A)= si nh*(A + „) ' (2 ' 6) 

Each solution R(X) of the Yang-Baxter equation fixes the structure constants of a Yang- 
Baxter algebra 

R 12 (X - yu)Ti(A)T 2 (/i) = T 2 (/i)T 1 (A)i? 12 (A - M ) (2.7) 

with generators T^(A), a,/3 = 1,2; where T\(A) = T(A) <g> /, T 2 (A) = I ® T(A) are the 
embeddings of the monodromy matrix T(A). 

Sklyanin's construction of open spin chains is based on the representations of two algebras 
T^ + ) and T^~> defined by the relations 

i? 12 (A - fx)T}-\\)R 12 (\ + /i)r 2 H (//) = Tt\v)Ri2{X + ^)T 1 ( - ) (A) J R 12 (A - //) (2.8) 



i? 12 (-A + / x)T 1 (+) ' 1 (A) J R 12 (-A -n- 2r ] )T 2 (+)t2 ( f i) = 

= T 2 Mt2 ^)R u (-\- f i-2r ] )T} +)h (X)R 12 (-\ + f i) . (2.9) 

We shall call and T* - ) right and left reflection algebras respectively. The transfer matrix 

t(X) =trT (+) (A)T ( - } (A) (2.10) 

as a trace in auxiliary space is the central object under consideration because it generates 
with [t( A), £(//)] = a commuting family of operators. 

The explicit construction of integrable open boundary conditions for models arising from 
the Yang-Baxter algebra starts with the 2x2 matrix 



K(\ n - 1 ( sinh(A + ne sinh(2A) \ 
1 ' U ~ sinh £ cosh A \Ke~ e sinh(2A) -sinh(A-0j 



, , 2/ve si nl ! A 2/vc ■ sinh A 

= I 2 + a tanh A coth£ H — — a + H — a 

sinh 4 sinh 4 

originally found by de Vega et al. [9]. It constitutes the known c-number representations 
K^- + \X) = |K(A + 77, £ + ) and K^~'(X) = K(X,^~) of the reflection algebras with the obvious 
properties 

trtf(A,0 = 2, K ( -\0)=I, trif (+) (0) = 1. (2.12) 

The Hamiltonian fll.ip involving two neighbouring sites for the interaction is connected, 
up to a factor, to the first derivative of t(X) by looking at the expansion t(X) = 1 + 2A7-^ + . . . 
around the point A = 0. Considering local L-matrices building up the two representations 
TW(A) = L L (X) ■ ■■L M+1 {\) and TH(A) = L M (X) ■ ■ ■ L X {X) of then by construction 



T { -\X) = T ( - ) (A)iT ( - ) (A)T ( - ) - 1 (-A) 
r (+)*(A) = TW(A)irW(A)(T( + )- 1 )'(-A) 



n t/ , (2-13) 
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are representations of the reflection algebras such that the normalized transfer matrix 

t(X) = K {+) (X)T(X)K { -\X)T- 1 (-X) , t(0) = l (2.14) 

is independent of the factorization of T(A) = T( + )(A)T(~)(A). Thus we are free to choose 

T (+) (A) = K {+) (X) , T { -\X) =T(X)K ( - ) (X)T-\-X) . (2.15) 

In order to gain more symmetric arguments and to avoid inconvenient scalar functions 
after applying the inversion formula 

r ' w ° ( W '-^ ^ (A - ,1 °' (2 ' 16) 

it is instructive to define the new object U(X + rj/2) = T^~\X) (d q T)(—X — 7]/2) consisting of 
U(X) =T{X-7]/2)K i -\X-7 1 /2)a y T t (-X-r ] /2)(x y . (2.17) 
It is still a representation of the left reflection algebra with a 2 x 2 matrix in auxiliary space, 

U(X) 



M(A) B(X)\ ( . 

lc(A) V(X) ■ (2 ' 18) 



The quantum determinant (d q T)(X) is the central element (Casimir) of the Yang-Baxter 
algebra (12.71) . With the one- dimensional projector P{ 2 onto the antisymmetric (singlet) state 
in the tensor product V <S> V of auxiliary spaces the definition reads 

(d q T)(X) = tr 12 P 1 ^T 1 (A - r]/2)T 2 (X + rj/2) 

= A(X + r ] /2)D(X-r ] /2)-B(X + i 1 /2)C(X-r ] /2) . l " ' 

Here, the trace tri2 is to be taken in both auxiliary spaces 1 and 2 of the tensor product 
V <8> V and the monodromy matrix T enters with the usual representation 



(A(X) B(X) 
1 ' \C{X) D(X) 



(2.20) 



3 Quantum Determinants 



Like the Casimir for the Yang-Baxter algebra there exists a similar object for the left reflection 
algebra. It is defined as 

(A q U)(X) = tr 12 P^U X {X - r?/2)i? 12 (2A - V )U 2 (X + V /2) . (3.1) 

To express (A q U)(X) in terms of the generators ^4(A), B(X), C(A) and V(X) it is instructive 
to use the combinations 

V(X) = sinh(2A)X>(A) - sinh rj A{X) , C(A) = sinh(2A + rj)C{X) (3.2) 
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borrowed from the algebraic Bethe ansatz. Then the suggestive form of the quantum deter- 
minant reads 



{A q U){X) = A{X + r ] /2)V{X-r ] /2)-B{X + r ] /2)C{X-r ] /2) . (3.3) 

In case of the c-number representation K(X — 77/2, £) for U(X) connected to the left 
reflection algebra the relation 

/a T ,w, / ^ sinh(2A — 2n) cosh A TW , ^ , 

(A q K) A - 77/2, £ = 1 hn " s det if A, £) 3.4 

cosh(A — 7/) 

holds. Note that this relation is only valid for the shifted argument A — 77/2 because the 
arising expressions are no longer of difference form. 

Appropriately transforming the boundary parameters £ and k according to the mapping 
[14] " 

sinh £ cosh £ 

sinh a cosh p 1 = , cosh a sinh f3 = (3.5) 

2k 2k 

the determinant det K (A, £) factorizes and its quantum version decomposes to product form 

(A q K)(X - 77/2, Q = - sinh(2A - 2 V ) = ^ ^ + a) "f ( ^ . 

smh a cosh p cosh(A — 77) sinh a cosh p cosh A 

(3.6) 

It is obvious that in the parametrization (13.51) the model is invariant under the simultaneous 
transformations a —>■ —a and /? — > i7r — /3. 

As the quantum determinant respects co-multiplication, applying it to the full represen- 
tation (I2.17P of the left reflection algebra with monodromy matrices T yields 

(A q U)(X) = (d q T)(X - 77/2) (A q K)(X - 77/2, £) (d q T)(-X - 77/2) . (3.7) 

Example. Consider the inhomogeneous periodic chain with inhomogeneities Sj G C at each 
lattice site j = 1 . . . L. Then the quantum determinant of a fundamental L-matrix Lj(X) = 
R(X — Sj) takes the scalar value 

{d q Lj)(X) = sinh(A - Sj + 377/2) sinh(A - Sj - 77/2) (3.8) 

for a chain of L local 



yielding (d q T)(X — 77/2) = nf=i srn h(A — Sj + 77) sinh(A — Sj — 77) 
spins 1/2. 
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4 Functional Bethe Ansatz 

The entries of the monodromy matrix U(X) can be expressed in terms of the operators A(X), 
B(X), C(A) and D(X) of the periodic monodromy matrix T(A), i.e. 

v ; cosh(A - r//2)sinh£- sinh(2A) v u ) \ u ) 

si „h ( A + ,/2 - r ) -h(2A - ,) _ ^ ( _ A _ 



cosh(A — ?y/2) sinh£~ sinh(2A) 

«Te sinh(2A — ri) . . „ . , 

+ — rr^ rr^ A - n/2) A(—X - n/2) 

sinh£-cosh(A-r7/2) v " ' v /7 ; 

_^^A-^) 

sinh£- cosh(A - r//2) v ;/ ; v ;/ ; v ; 

Obviously there is no easy pseudo vacuum |0) for an algebraic Bethe ansatz to work (23(A) |0) = 
0) for non-diagonal boundaries (k~ ^ 0) in K^~\ Instead we can apply Sklyanin's functional 
Bethe ansatz [17]. 

Henceforth let us restrict to the rational case. Then we can choose the boundary matrices 
diagonal in favour of a twisted monodromy matrix T(A) as demonstrated in the following 

Example. Applying the rational limi^\ to the parametrization of the boundary matrix K(X,£) 
yields the similarity transformation independent of the spectral parameter 



K(X,0 = Ms( a « X ^(MSy 1 (4.2) 



with the diagonal matrix to be a solution to the reflection algebra and the 2x2 number 
matrices 

' e +6/2 o \ l / e /3/2 _ e ~/3/2> 



Thus in the transfer matrix t(X) = tr K^(X)T(X)K^(X)T- 1 (-X) we are free to consider a 
diagonal outer boundary matrix together with the c-number twist (M^ S^ + ')~ 1 T(X) of 
T(X)E 



Operator- Valued Zeroes 

The main goal of the functional Bethe ansatz is to treat the spectral problem of the transfer 
matrix in a representation space of symmetric functions manipulated by some shift operators, 
which descent from operator-valued zeroes of the £>-operator. Starting from B(X) in the 

1 'Rational limit' means to scale the arguments by a small parameter e — > and calling r/ — > sic. Here 
only a is rescaled whereas [3 remains unchanged. This immediately follows from (|3 - 5[) . 

2 The gl(2) symmetry of the rational model allows to remove the twist of the monodromy matrix in favour 
of a twisted boundary matrix = (M<+) S^y 1 AfW 5^, see Ref. [12]. 



s 



rational limit reading 

B(X) = 2A " 1C 



A ! 4 1C/ -B(-X - ic/2)A(X - ic/2) 



2A 

-A + £- - ic/2 



B(X - ic/2)A{-X - ic/2) 



-2X 

- K-e ~A(X - ic/2)A{-X - ic/2) 
+ K~e' e 'B(X - \c/2)B{-\ - ic/2) 



(4.4) 



we observe the expression in the square brackets to be symmetric with respect to A — > —X 
and having no pole at A = 0. As the operators A, B, C and D are polynomials^ of degree L 
with the known asymptotics 

V2 cosh (3+ V2 cosh (3+ 

V2 cosh V2 cosh^+ 

we can factorize the square brackets in terms of A 2 . The asymptotic pre-factors arise from 
the twist S~ l M _1 of the periodic monodromy matrix T(A). 

Thus B(X) is polynomial with a simple zero at A = ic/2 and operator- valued coefficients 
assembling 

As [23(A), $(//)] = 0, according to the reflection algebra, we can deduce [x 2 ,^ 2 .] = for all 
j, k = 1 . . . L with the spectrum shown in the next 

Example. Setting L = 1 in the explicit expression (14. 6 j) /or a spin 1/2 representation with 
inhomogeneity s\ yields for the argument A = the form 

^( (Sl + ,C/2)2 (.-.^) ^ 

o/ i/ie operator-valued zero x\ reading in a diagonalized form. This fixes discrete sets Aj = 
{sj — ic/2,Sj + ic/2} representing the spectra of the coordinates Xj except for a global sign 
and resembling X L = Ax x . . . x Al- 

For the forthcoming relations, to work with the simple zeroes Xj instead of x 2 , we refer 
to the supplemental 

Remark. Let Xj be the operator-valued zeroes satisfying _\=±% \B(X) = 0. Then all x 2 can 
be simultaneously diagonalized such that for all j,k = 1 . . .L 

[ Xj ,x k ]=0 , x 3 =^ + iG/2 s ._ ic/2 ) ■ (4.8) 



3 From here the e^-dependence of all periodic chain operators is suppressed. 
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Conjugated Momenta 

With the operators Xj the next problem is to calculate the expression for the transfer matrix 
in the x-representation. For this let us introduce first the 'conjugated momenta' to the 
'coordinates' Xj. 

Considering A(X) and T>(\) as polynomials and inserting the operator valued zeroes of 
£>(A) by 'substitution from the left' yields the new operators 



V 



£>(A)=$>^X+ 



(4.9) 



Here A p and V p denote operator-valued expansion coefficients. The commutation relations 



with the coordinates j £1X6 summarized in 

Theorem 1. Letxj be the operator-valued zeroes of B(X) and their conjugated momenta 
related by the reflection algebra. Then 

Xfx k = {x k ±ic8 jk )Xf . (4.10) 

Proof. Consider the commutation relation 

\ sinh(A + fi - rj) sinh(A - \i - rf) 

A(X)B(fi) = —T7T \ ■ i h ^ — -B(/j,)A(X) 

v ' smh(A + fj) smh(A - //) 

sinh?7sinh(2/i — ri) ^ 

+ . - /n ' , . \ f. 'ArB(\)A(n) 4.11 

sinh(2/x) smh(A + fi) v ' K ' K ' 



SinhT] -B(X)V(») 



sinh(2/x) sinh(A + fx) 



of A(X) and £>(/i) and multiply it by sinh(A + fj) sinh(A — //). Then take its rational limit 
and insert the coordinates Xj by 'substitution from the left'. The expression reduces to 



x 



: j - X)(xj + X)X-B(X) — (x + X- ic)(x - A - ic)B(X)X~ . (4.12) 



Replacing B(X) by its factorized form, cancelling the constant asymptotics and multiplying 
by the inverse (x? — A 2 ) -1 from the left the commutation relation 

X i[fl( X2 -^} = [A 2 -(%-ic) 2 ][n(A 2 -^)]xr (4.13) 
i=i i=i 

remains. Implying all expressions to be symmetric in Xj to act on we arrive at the desired 
relation. Analogously the elementary commutation of with the coordinates arises from 

the commutation of B and V. □ 

The next natural step would be establishing the commutation relation between two X^s. 
However, it cannot be done directly because Xp exceed per definition the representation 
space. 



10 



Representation space 

Following Sklyanin's approach the square brackets in operator (14.61) can be expanded into 
\ 2L — bi\ 2(yL ^ ± . . . + bi with commuting operators bj thus sharing a common system of 
eigenfunctions f a , 

%f« = Vjf a , a = 1...2 L (4.14) 

where 2 L represents spin- 1/2. To every point h a = (b", . . . , 6°) G B L C C L there corresponds 
only one eigenfunction f a and the representation space W (e.g. for the XXX chain we have 
W = (C 2 )® L ) is isomorphic to the space FunB L . 

Example. A possible realization of the eigenfunctions f a is 

$ j f a )Qf)=qf a {bP) (4.15) 

where bj act as multiplication operators. Let {t a G (C 2 )® L | (t a )P = / a (b^)} be a basis ofW 
then, with the constraint / a (b*) = 8@, it is indeed orthonormal and complete. 

Since b n are the symmetric polynomials of the roots x 2 we are led to consider the mapping 

6:C L ^C L , x^b (4.16) 

given by the formula b n (x) = s n (x). The s n (x) are the elementary symmetric polynomials 
of degree n — 1 . . . L of c-number variables 

si (x) = x\ + x\ + . . . + x\ 

'■■ (4-17) 

/ \ 2 2 2 

S M X J — X^X 2 • • • x L 

The diagram 

{0,1}CC 




of the combined mapping / o 9 reveals the isomorphism between Fun B L = W and the space 
of symmetric functions SymFunX L . The set {0, 1} is the range of / in the example above. 
Thus the operator roots x 2 can be thought of multiplication operators 

&j9(Vi, ■■■,Vl) = Vjgivi, ■■■iVl) (4.18) 

in an extended representation space FunX L = W which is a non physical one. Recall all the 
results should only use the original space SymFun X L = W, as the operators A, B, C and V 
map SymFun X L — > SymFunX L or more sloppy W W . So for plainness we will use in the 
following the terms W and W for the representation spaces instead. 

For the action of X- on a function s G SymFun X L we need first to extend the operators 

from W to W by the constant function 

w(x) = 1 for all x G X L . (4.19) 
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Obviously u is symmetric and thus belongs to the representation space W. Now define the 
action of Xj on u by 

[Xfu,)(*) = &f(x) . (4.20) 

Then the functions A - (x) uniquely determine the action of Xj on any vector s which is 
identified due to the isomorphism with some symmetric function s(xi, . . . ,xl) = (su)(x) 
created from the cyclic vector uj by the operator 's = s(xi, . . . , xl). Thus 

(^)(x) = (Xf^)(x) = S (£±x)(X±u;)(x) = x)A±(x) . (4.21) 

Here we introduced the shift operators 

Ef : C L -> C L : (xi, . . . , Xj , . . . , x L ) i-> (a*, . . . , x,- ± ic, . . . , x L ) (4.22) 

acting on some L-tuple of c-numbers. In the extended representation space W of not neces- 
sarily symmetric functions the action of X- then reads 

Xf = AfE± (4.23) 

with A^ = A^(x). By the operator relation (14.231) we can now calculate the commutations 
of the momenta in 

Theorem 2. Let Xf be the conjugated momenta related to the coordinates Xj by the reflection 
algebra. Then 

\Xf,X£]=0 for allj,k = l...L , 

L 3 J (4 24) 

[Xf, Xj;] = for all j,k = l...L but j ^ k . 

Proof. Let us start with X~ where the first assertion is obvious for j = k. Then it is enough 
to consider the cases j = 1, k = 2. Taking the rational limit of 

[A(X),A(ji)] = n) ^ C(A) " B{X) C(/i) ] (425) 

and inserting A = x\ and \i = Xi by 'substitution from the left' the RHS turns into zero and 
for the LHS we get 

\=x\.p=X2 A(\^A(fl) ^ X\ X2A m A n ^ ^ X>2 X] AmA n 

m,n m,n ^ ^ 

= ^ x 2 X± A n = Xi ^ ^ x 2 ^4 n = -X^ ^2 

n n 

In the same way starting from A(fJ,)A(X) one obtains X 2 ~X 1 ~ and the assertion is proven. 
The commutation of X + 's and mixed commutators excluding the cases j = k can be treated 
analogously by considering 

V(\),V{n) = - sinh(2A + rj) sinh(2/j + 77) [.4(A), .4(//)] 

«(A). ,%)] = g^±g||g|^ A{p)] < 4 ' 27 > 

in the rational limit. □ 
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The remaining commutation involving the quantum determinant A q is summarized in 

Theorem 3. Let Xj and Xf be the coordinates and conjugated momenta related by the 
reflection algebra and A q (X) is the quantum determinant. Then 

XfXj = A q (x j ±ic/2)forallj,k = l...L . (4.28) 

Proof. Substituting the operator-valued zeroes Xj into the quantum determinant (j3.7p one 
obtains 

A q ( Xj - ic/2) = ~ ic T A rnV n = ~ ic) n {xJA m )v n 

m,n m,n 

= J2@J - K) n XjV n = XJ ( 4 ' 29 ) 

n n 
= X7X+ 

and analogously A q (xj + ic/2) = X^XJ exerting the reflection algebra. □ 

Remark. The remaining zero ic/2 ofB(X) is an exception and renders the operators A(ic/2) = 
d q (—ic/2) and V(ic/2) = to be constant yielding A q (ic) = 0. 

Representation of A ± 

Applying XfXj = A^E^AJEJ to an arbitrary function g G W induces the sequence 

{XfXjg)^) = {AfEfAjEjg){*) = A±(x)(£±Aj£f ^(x) 

= Af(x)(AjEfg)(Efx) = Af(x)AJ(Efx)(E?g)(Efx) 

= A±(x)AJ(£±x)«7(x) (430) 

= A q ( Xj ± ic/%(x) 

relating the representations A^ to the quantum determinant A q . In the case of a finite 
dimensional representation of the generators {xj, Xj t }j !/ =1 such that the spectrum X L shows 
no multiple points the problem of constructing such a representation is equivalent to that of 
determining the functions {A^}^ =1 on X L satisfying 

A±(x)A±(£±x) = A±(x)A± (£±x) for all n, m 

A+(x)A;(E+x) = A;(x)A+(^x) for all n,m but n ? m (4.31) 
A, fa ± ic/2) = Af (x) A J (£±x) for all j 

arising from theorems [2] and [31 The above relations are not defined when the shifts E^ move 
the point x out of X L = Ai x . . . x A^. This means {A^}^ =1 have to vanish on the boundary 

dXf = {x G X L |£±x G C L \X L } (4.32) 

of the set X L . For the open XXX chain with Aj = {sj — ic/2, Sj + ic/2} this is clear from 
the following 
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Example. The vanishing o/A^(x) on the boundary <9X^ can be directly seen from the explicit 
factorization of A q (X) = A~(A + r]/2)A + (X — r\/2) into 



A~( A j = A t lC H^ a \ - si + ic/2) (A + Sl + ic/2 



-1 L cr 



A+(A) = -(2A - ic)e X+ / C l* a \f[(X - «, - ic/2)(A + Sl - ic/2 



(4.33) 



i=i 



considered in the rational limit indicated by e — > 0. 



Spectral Analysis 

Now let us return to the original problem, the spectral analysis of the new transfer matrix 
r(A) = trK(\ + rj/2,£ + )U(\)/2 in the rational limit 

(A + ic/2)(A + e + -ic/2) l A-e + + ic/2 ~ M , 
H A ) = ^ ^( A ) " 4^|+ V W 



+ 



2A£+ 
(A + ic/2)«" 



e^C(A)+e- 9+ S(A) 



(4.34) 



and mind the scaling factor e — * 0. To plug in the zeroes Xj by 'substitution from the 
left' we have to get rid of C(A) by diagonalizing K^ + \ Thus only the first line remains. The 
diagonalization does not change the quantum determinant A g (A) and the eigenvalue problem 
r(X)(f = A(\)ip can be solved by 'substitution from the left' reading 



\=Xj 



|r(A) 



(xj + ic/2) (xj + q + -ic/2) ^_ _ l %-q+ + i£/2 x . 



2xjct~ 



e 4xja+ 



(4.35) 



With this representation at hand one observes 'separation of variables' suggesting the product 
ansatz 



l[Q(xi) 



(4.36) 



z=i 



for the eigenfunction (p e SymFunX L = symmetric in its arguments To explicitly 
apply the operator-valued expression (I4.35P one should clarify its behaviour by the following 

Lemma. The action of the combined expression XjXj by 'substitution from the left' onto a 
symmetric function s = s(xi, . . . , xl) is given by 

XjXfsix) = (%X±s)(x) = Xj (Xfs){x) = Xj Af{x)s{Efx) . 



(4.37) 



Then applying (14.351) to ip only the jth argument is affected such that the problem 
separates and 



a/ \ / x (xj + ic/2) (xj + a + - ic/2) 



2 X j Gi ~^ 



1 Xj — a + + ic/2 



(4.38) 



Axj(x + 



A + (x j )Q(x j + ic) 
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holds. Here we used (I4.37P with the allowed arguments Xj G Aj = {sj — ic/2, Sj + ic/2} on 
the grid entering A - (x) = A ± (xj). 



Remark. Along with the normalization (d q T)(—X), c.f. (13.81) . and a shift in the arguments 
the original transfer matrix t(X) is related to r(A) via r(A) = (d q T)(— A) t(X — T}/2). 

5 TQ-Equations 

The eigenvalue problem as formulated in Eq. (I4.38P reduces to a system of homogeneous 
linear equations due to the fact that A ± (xf) = at the points xf = Sj ± ic/2: 



(xf + ic/2)(xf + a+ - ic/2) 
A(x+)Q(x+) = ^ ^A-(x+)Q(xJ 



A(xj)Q(xj) = — 3 AxTa+ ' ^(xj)Q(xj 



(5.1; 



For pairwise different inhomogeneities, Sj ^ for j ^ k, these linear equations allow for a 
non-trivial solution provided that the following functional equation for the eigenvalues A are 
satisfied^ 



Afe + ic/2)Afe - ic/2) = -a±^ ( j;: l )a+ £1 ° )a+ A a fe) . J-I...L (5.2) 



Using the known asymptotic form of the even polynomial A(A) = A(— A) we are led to the 
ansatz 

.... (— l) L sinh/3 + sinh/3~ + cosh(# + — 9~) 2J , 2 ^ 2r .,, 2 „. 

A A = ^— J- — i '-X 2L+2 + a 2L X 2L + a 2L _ 2 X 2L - 2 + ... + a . 5.3 

a + a~ cosh p + cosh p 

The (I/ + 1) unknown coefficients aj are determined by Eqs. (15.21) and the constraint A(ic/2) = 
d q (—ic/2) with the quantum determinant d q (X) of the periodic chain. Thus the solution of 
the spectral problem amounts to finding the common roots {ct>2j}f=oi v = 1 . . . 2 L , of these 
polynomial equations. This task is of the same complexity as finding the eigenvalues of the 
spin chain Hamiltonian directly and therefore this approach is limited to small system sizes 
where we have checked numerically that it does indeed yield the complete spectrum. 

To compute the eigenvalue of the transfer matrix or the spin chain Hamiltonian in the 
thermodynamic limit L — ► oo the functional equations introduced above need to be analyzed 
beyond the set X L using explicitly the analytic properties of the functions therein. 

Remark. Note that the functional equations $5.2\) hold at the discrete points Sj only: treating 
the Sj as a continuous variable and applying standard Fourier techniques one can compute 



4 If n of the inhomogeneities coincide the (n — 1) derivatives of this equation at this value of Sj have to be 
taken into account in addition. 



15 



lnA(A) and thereby the corresponding eigenvalue of the spin chain Hamiltonian U.l\) , i.e. 
TC = ic (d\hxr)(ic/2) = ict'(O) in the homogeneous limit. For \a \ > c/2 one obtains 

ic^^(ic/2) =^(\a + \/2c) - ij(\a + \/2c + 1/2) + c/\a + \ 

dX 

+ ip{\a-\/2c) - ip{\a-\/2c + 1/2) + c/\a~\ ( 5 - 4 ) 
+ 7r-21n2-l + (2-41n2)L 

which is for imaginary a the known energy eigenvalue of the XXX spin chain with diagonal 
boundary fields [10] (ip is the dig amma function). However, the non-diagonal contributions 
and corrections of the order 1/L are not included. This is a consequence of neglecting the 
corrections to Eqs. ( E3| ) away from the points Sj . 

Instead we go back one step and consider the Eqs. (14.381) for general arguments Xj — > A. 
Formally, this is a second order difference equation reminiscent of Baxter's TQ-equation 
[2]. The analysis above leading to Eqs. (15.21) has been based on the singular points of the 
TQ-equation at the boundaries <9X L , i.e. points where one of the coefficients A vanishes. 
Away from these points there exist two independent solutions to (I4.38P and one needs some 
information on the properties of the unknown functions Q(X) in this formulation of the 
spectral problem which has to be solved for polynomial eigenvalues A(A) of the transfer 
matrix. 

In cases where a pseudo vacuum exists and the algebraic Bethe ansatz is applicable 
to the solution of the problem the Q-functions are known to be symmetric polynomials 
Q(X) = n^=i(^ ~~ v e)(X + vi) with roots ve satisfying Bethe ansatz equations. Note that only 
in these cases the constant function u — 1 introduced in the construction of the representation 
of the Yang-Baxter algebra on the space W can be identified with the pseudo vacuum |0). 

In general, the TQ equation can be rewritten as a recursion relation 

Q(X + ic) = o(A)Q(A) + 6(A)Q(A - ic) (5.5) 
for the function Q(X) or equivalently, with the auxiliary function P(X + ic) = Q(X), 



Q(X + ic)\ _ fa(X) b(X)\ (Q(X)\ 
P(X + ic)) \1 ) \P{X)) 



(5.6) 



The coefficients a(A) and 6(A) are obtained from the TQ-equation (14.381) and show constant 
asymptotics for large values of their arguments 

A(A) AeXa+ o sinh/3 + sinh/?" + cosh(#+ - 0~) 

a \ X ) = — ~a~Z7\T \ 1 , ■ /o ~ 2 



A+(A) A-a+ + ic/2 cosh (3+ cosh (3~ 

b(X) = A-(A)2(A + ic/2)(A + a+-ic/2) g _ ( ' j 

1 ' A+(A) A-a+ + ic/2 ~ ' 

This allows to solve the recursion relations in the asymptotic regime |A| 3> 1 yielding 

\n \n \n— 1 \n— 1 

Q(A + inc) = 7 , 2 Q(X + ic) - A X A 2 ^- ^_Q(A) . (5.8) 

Ai — X2 Ai — A2 
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Here n is an integer and Ai^ = e ±( ^ are the eigenvalues of the asymptotical matrix of coef- 
ficients in (15.61) . cosh0 = (sinh f3 + sinh f3~ + cosh(0 + — 6~ )) /(cosh (3 + cosh (3~). Ordering the 
eigenvalues as |Ai| > | A2 1 we obtain for fixed A and large n the leading term Q(X + inc) ~ A" 
suggesting the following ansatz for the asymptotic form 

,A0 



Q(X) ~ exp (t-) x polynomial in A 



IC 



(5.9) 



Here the polynomial form of the subleading part assures that the eigenvalue A(A) of the 
transfer matrix remains polynomial. Note, that since A(A) is an even function of its argument 
there exists always a second solution Q(— A) to the TQ equation which decays exponentially 
for A — ► 00. 

Only in two cases, namely = and ire or, equivalently, 

cosh(#+ - Q-) = ± cosh(/3+ t P~) , (5.10) 
the exponential factor disappears and the TQ equation can be solved by an even polynomial: 
in the first case Eq. (I5.9P implies that Q(\) = n*=i +) (^ — v e){^ + v i)- For <p — ' m , the 
exponential factors can be removed by the transformation Q(X) = exp(iA7r/ic)<5'(^) resulting 
in a TQ equation for Q'\ 

(x J +ic/2)(x j + a + -ic/2) A _ / 



A(xj)Q'(xj) = - 



2^2/ j Ot 



+ 



a 



+ ic/2 



4:Xja~ 



-A (xj)Q'(xj — ic) 
^ + {x j )Q'{x j +\c) . 



(5.11) 



Again, it follows from the asymptotic analysis that this equation allows for a polynomial 
solution Q'(X) = n^=i \^ ~ v t)(^ + v e) whose existence has been verified by numerical 
analysis for small systems. 

In both cases the spectrum is determined by the roots of these polynomials. To guarantee 
analyticity of the transfer matrix eigenvalues A(A) the Vj, j = 1 . . . have to satisfy the 

Bethe ansatz equations 

L 



v j + a — ic/2 vj + 



ic/2 



Vj — a~ + ic/2 Vj — a + + ic/2 



1=1 ° 



si + ic/2 Vj + si + ic/2 



si — ic/2 Vj + si — ic/2 

r M(±) 

n 



v k + ic Vj + v k + ic 



fc=i 



v i — Vk — IC Vj + V k — ic 



(5.12) 



Note that Eq. (15. 10|) is equivalent to the constraint that the boundary matrices can 
be simultaneously diagonalized or brought to triangular form. In this case Eqs. (15.121) can 
be obtained by means of the algebraic Bethe ansatz [12] or in the rational limit from the 
TQ-equation approach for the open XXZ chain [19]. In this trigonometric case the complete 
set of eigenvalues is obtained from two sets of Bethe equations which both reduce to (I5.12p . 
This is due to the invariance of the model under the change of parameters a — > —a and 
(3 —>■ iir — j3 which maps <p = ire, see Eq. (13.51) . As another difference to the situation 
in the XXZ model the number of Bethe roots is not restricted by the constraint on the 
boundary fields and we have to consider solutions of the TQ-equations 'beyond the equator', 
M(±) > L/2. 
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6 Summary 

In this paper we have extended Sklyanin's functional Bethe ansatz method to the open XXX 
chain with non-diagonal boundary fields. Within this framework we have derived a single TQ- 
equation (14.381) which determines the spectrum of this model for any values of the boundary 
parameters. This supports the approach of Yang et al. [19] who obtain a TQ-equation for 
the XXZ model assuming the existence of a certain limit in the auxiliary space. 

The TQ-equation allows for a solution in terms of polynomials for the function Q provided 
that a constraint between the left and right boundary field is satisfied. In this case the solution 
is parametrized by the roots of one set of Bethe ansatz equations. In our numerical study of 
small systems we obtain the complete spectrum from these equations, when solutions 'beyond 
the equator', i.e. polynomial Q with degree > L/2, are taken into account. This is in contrast 
to the trigonometric case where two types of Bethe equations are required and the degrees 
of the corresponding Q-functions are fixed by the constraint. We suppose that this feature 
of the XXX case is a consequence of the rational limit. 

If the constraint between the boundary fields is missing only the asymptotic (exponential) 
behaviour of the Q-functions is obtained from the TQ equation, the subleading terms have to 
be chosen such that the eigenvalues of the transfer matrix remain polynomial. In Sklyanin's 
approach the Q-functions contain all the information on the eigenstates of the model. There- 
fore, their determination for generic boundary parameters is necessary to tackle the problem 
of computing norms and scalar products within this approach. To make progress in this 
direction it should be useful to investigate how the recent construction of Galleas [11] con- 
nects to the TQ-equation (I4.38p . His solution of the spectral problem for XXZ chains with 
non-diagonal boundaries is given in terms of the zeroes of the transfer matrix eigenvalues 
and two complementary sets of numbers which parametrize matrix elements of certain ele- 
ments of the Yang-Baxter algebra. They satisfy equations reminiscent of the nested Bethe 
ansatz used to solve models of higher-rank symmetries. Further studies are necessary to see 
whether this parametrization of the spectrum in terms of O(L) complex numbers can be 
used to obtain a closed expression for generic (non-polynomial) Q-functions. This would be 
of great importance for the applicability of the TQ-equation to solve the spectral problem of 
integrable quantum chains. 
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